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After a brief introduction to Heavy Quark Effective Theory, we discuss a representation 
in HQET and methods of calculation of some kinds of HQET diagrams up to three loops. 



1 Introduction 

Effective field theories are very useful for describing physics at low energies ^ M, or large 
distances S> 1/M, where M is a high energy scale where some new particles or interactions 
become important. Effective Lagrangians are constructed as series in 1/M. Coefficients in them 
are obtained by matching scattering amplitudes in the full theory and in the effective one up 
to some order in 1/M. These matching coefficients are the only quantities which depend on 
M. All calculations inside the effective theory involve only characteristic energy scales (they 
are ^ M) of processes under consideration. The case when there is one such scale is especially 
simple. We can choose the renormalization scale n of order of this characteristic energy scale. 
Then there will be no large logarithms in perturbative series, and truncating such series will 
produce small errors. If we try to consider the same process in the full theory, there is a second 
scale M, and no choice of fj, allows us to get rid of large logarithms. Also, each extra scale in 
Feynman diagrams (with loops) makes their calculation much more difficult technically. 

Heavy Quark Effective Theory (HQET) is an effective low-energy field theory for some 
problems in QCD. In many widely-known effective field theories (Heisenberg-Euler theory of 
low-energy photon interactions, Fermi 4-fermion theory of weak interactions at low energies) 
the heavy particle (electron or W in these examples) does not appear. In HQET, the heavy 
quark appears in initial and final states, but is always nearly on-shell and non-relativistic (in 
some reference frame). 

HQET is discussed in textbooks [U [2] in detail. Here we shall concentrate on methods and 
results of calculations of multiloop Feynman diagrams in HQET, see e.g. [3]. Various methods 
of multiloop calculations are presented in the excellent book |1] in great detail, most of these 
methods are used in HQET. 



2 Heavy Quark Effective Theory 
2.1 Lagrangian and Feynman rules 

We are going to consider a class of QCD problems involving a single heavy quark with (on-shell) 
mass 771 3> Aqcd- Namely, we require that there exists a reference frame where it stays nearly 
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at rest all the time. In other words, there exists a 4-velocity v {v^ = 1) such that 

p = mv + p , (2-1) 

and the characteristic residual momentum <^ m. Light quarks and gluons also have char- 
acteristic momenta p^ <C m. Such problems can be described, instead of QCD, by a simpler 
effective field theory called HQET. Its Lagrangian is a series in 1/m. At the leading order, 

L = QJv ■ DQ^ + O + (light fields) . (2.2) 

The HQET heavy-quark field satisfies fQ^, = Q„. All light fields are described as in QCD. In 
the V rest frame, 

L = Q+iDoQ + O (^-^^ + (light fields) , (2.3) 

where Q is a 2-component spinor. 

The mass shell of the heavy quark, i.e., the dependence of its residual energy po on its 
momentum p, is 

Po^Pa-m . 

2m 

At the leading order in 1/m, it becomes po = 0. This is exactly what follows from the La- 
grangian (j2.3p . 

The HQET Lagrangian p.2p is not Lorentz- invariant, because it contains a fixed vector v. 
However, v is not uniquely defined. It can be changed by ~ p/m (sec (j2.1[) '). Lagrangians with 
such different choices of v must produce identical physical predictions. This requirement is 
called reparametrization invariance, and it restricts 1 /to" corrections in the Lagrangian. 

The heavy-quark chromomagnetic moment is, by dimensionality, ~ l/m. Therefore, at the 
leading order the heavy-quark spin does not interact with the gluon field. We may rotate the 
spin at will without changing physics — heavy-quark spin symmetry. In particular, B and 
B* are degenerate and have identical properties, because they can be transformed into each 
other by rotating the b spin. We can even change the magnitude of the heavy-quark spin (e.g., 
to switch it off) without changing physics — this supersymmetry group is called superfiavour 
symmetry. 

It is difficult to simulate a QCD heavy quark on the lattice because the lattice spacing 
a must be much less than the minimum characteristic distance of the problem, 1/m. The 
HQET Lagrangian does not contain to, and the only applicability condition of its discretization 
is a <C l/p- When we investigate the structure of heavy-light hadrons, p ~ Aqcd, and the 
condition a <^ 1/Aqcd is the same as for light hadrons. 

The Lagrangian (|2.2p gives the Feynman rules 

a fi 

^ = z^ - \ > ^ > =»g^°^r (2.4) 



p 2 p ■ V + iO ' 

In the V rest frame, the propagator is (the unit 2x2 spin matrix assumed) 

1 



Po + «0 
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In the coordinate space, the heavy quark docs not move: 







-i9{xo)S{x). (2.6) 



HQET-quark loops vanish because the heavy quark propagates only forward in time. We can 
also see this in the momentum space: all poles of the propagators in such a loop are in the 
lower po half-plane, and closing the integration contour upwards, we get 0. 

These Feynman rules can be also obtained from QCD at to — > oo. The QCD massive-quark 
propagator gives the HQET one: 



O 



P 



mv + p p ym 

m + m^+p 1 + 1 . QfP_\ (27) 



TO 



(mv + p)^ — m'^ + iO 2 p ■ v + iO 
The QCD vertex, when sandwiched between two projectors, becomes the HQET one: 

i±l ,i±i^l±iy_±l^ (2.8) 

2 ' 2 2 2 ^ ' 

When there is an external leg near a vertex, there is no projector; but we can insert it, and the 
argument holds. 

We have thus proved that at the tree level any QCD diagram is equal to the corresponding 
HQET diagram up to 0{p/m) corrections. This is not true at loops, because loop momenta can 
be arbitrarily large. Renormalization properties of HQET (anomalous dimensions, etc.) differ 
from those in QCD. 

2.2 One-loop propagator diagram 

Let's calculate the simplest one- loop diagram (Fig. [T|) 

1 r df^k 



J ^_2(k+p)-v- iO] [-fc2 - iO] 

1 f dknd'^-^k 



n2 



•^/^ / [_2(fco+^)_^o]"^[-fc2 + fc2_,o]"^ 



^{-2ujf-'''~^''U{ni,n2). (2.9) 



It depends only on the residual energy lu = po, not p; the power of -~2uj is clear from dimensional 
counting. 

If > 0, real pair production is possible, and we arc on a cut. We shall consider the case 
uj < 0, when the integral is an analytic function of ui. We'll set ~2uj = 1. If ni is integer and 
n-i < 0, /(«!, ^2) = because this is a massless vacuum diagram. If 712 is integer and n2 < 0, 
I{ni, 712) = because the diagram contains an HQET loop. 

At w < 0, all poles in the kg plane are below the real axis at fco > and above the real axis 
at fco < 0, and we can rotate the integration contour counterclockwise without crossing poles 
(if > 0, we cross the pole at ko = —uj — iO). This Wick rotation 

ko = ikEo (2.10) 
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k + p 



Figure 1: One-loop propagator diagram 



ko 

'\k\+iO 



\ C 



uj-iO \k\-iO 



Figure 2: Wick rotation 



(Fig. [2|) brings us into Euclidean momentum space (fc^ 
In the Euclidean space, 



-4) 



I{ni,n2) 



Using the well-known formula 



1 



dk 



EO 



d'^-^k 



{k^ + kl 



EO) 



T{n) 



with d d — 1, m? — > A;|;q, n — > 712, we obtain 



(2.11) 



(2.12) 



/(71l,7l2) 



7rl/2r(7l2) 



(fc|;o)' dkEO 



l-2ik 



EO 



(2.13) 



The integrand is even in kso, and has cuts at k^Q < (it would be wrong to write it as 



^d-l-2n2^ 

the integral via the discontinuity at this cut: 



k^o^ ^"^). Deforming the integration contour around the upper cut (Fig. [3|), we can express 



r(n2 



J(ni,n2) = 2 \ \ 
7r^/^r(n2) 



■ cos 



00 ^d-l-2„2^^ 

(2fc + l)»i 



(2.14) 
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Figure 3: Integration contour 



This integral can be easily calculated in F functions: 



7(711,712) 



2^2712 —1^+1 

7ri/2 



■ cos 



Tjd - 27i2)F(ni + 2112 - d)T{n2 - ^) 
F(ni)F(n2) 



Using the well-known properties of the F function 

F(2a:) = 7r-i/222--ir(a;)F(x + i) , F(x)F(l - x) 
we can simplify this result: 

F(7ii + 2n2 -d)F(f -712) 



sm nx 



7(711,712) 



F(ni)F(7i2) 



(2.15) 



(2.16) 



(2.17) 



It is also easy to derive this result in coordinate space [3] . HQET propagators in momentum 
and coordinate space are related by 



duj 



it 

{-2uj - iO)" 2tt 2F(77.) V 2" 



,(iw-0)t 







it 



dt 



2iT{n) 



(-2cj-i0)" ' 



massless propagators — by 



i T{d/2-n) 



d/2-n 



(-p2 - iO)" (27r)'* (47r)'^/2 T{n) \-x^+iQ, 
4 V .... ... .w/oF(d/2-n) 1 



-x"^ + iO 



) r(n) 



(-p2 - i0)''/2-" 



(2.18) 
(2.19) 

(2.20) 
(2.21) 



Our diagram in coordinate space (Fig. 31 x = vt) is just the product of the heavy propaga- 
tor (Pl^ and the light one ^iW^ (where -x^j^ = -t^/A = {it/2f): 



1 1 F(d/2-n2) (it 



2 (47r)'^/2 F(ni)F(7i2) V 2 



ni+2n2 —d—1 



0{t). 
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ni 



Figure 4: One-loop propagator diagram in coordinate space 

The inverse Fourier transform (j2.19p gives our diagram (j2.9[) in momentum space 

^/K,n2)(-2a;)'^-— , 
where /(ni,rt2) is given by (|2.17p . 
2.3 Renormalization 

The Lagrangian contains bare fields and parameters: 

L = Q^^iv ■ DoQvo Do^ ^ - tgoA^^t" . (2.22) 
They are related to the renormalized ones by the renormalization constants: 

Q„o = ^q^'Q., ^o = ^y'A, ao = ZAa, go = Z'J^g, (2.23) 

where ao is the gauge-fixing parameter. The ghost field, the light-quark fields (and their masses) 
are renormalized as in QCD (not written here). Minimal renormalization constants have the 
structure 

^11 Os ^ / Z22 ^ Z2i\ /Os^ 2 
47r V ^2 £ ] V47r^ 



l + + +... (2.24) 



They don't contain and e" {n > 0) terms, only negative powers needed to remove divergences, 
and hence are called minimal. We have to define as to be exactly dimensionless. In the MS 
scheme as depends on the renormalization scale ^: 

d ^ ^2e^zMt^))e''^' . (2.25) 



(47r)'i/2 ^ 47r 

where je is the Euler constant. 

Let's calculate the HQET propagator with one-loop accuracy: 



+ 



iSiuj) = iSoiuj) + iSoiuj)i-i)T.{uj)iSoiuj) + ■■■ (2.26) 

where 

So{lo)^-. (2.27) 

Ul 
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The one-loop heavy-quark self-energy (Fig. [5|) is 



S(ti;) = iCf 



(Pk 



ko + cJ'^'"'^ k^ ^ P 



(2.28) 



where ^ = 1 — qq. In the numerator, we may replace (fc • v)'^ = (kg + w — lo)'^ — > w^, because if 
we cancel kg + u> in the denominator the integral vanishes. Using (|2.17p . we obtain 



I](w) = Cf 
= Cf 



2/(l,l) + |/(l,2) 



(47r)d/2 

ffg(_2c.)i-2^r(l + 2e)r(l-£) 



d-4 



d-3 



(2.29) 




Figure 5: One- loop heavy-quark self-energy 
The propagator expressed via renormalized quantities is 



ujS{uj) = 1 + Cf- 



47re 



(2.30) 



where 



L = log 



-2cj 



It should be equal S{uj) 
Therefore, 



ZgSrii^), where the renormalized propagator Srit-^) is finite at e ^ 0. 



1 



.CH3-a)^ 



(2.31) 



(it is also easy to write Sr{uj)). 

The HQET field does not renormalized in the Yennie gauge a = 3. This is exactly the 
reason why this gauge has been introduced in the theory of electrons interacting with soft 
photons (Bloch-Nordsieck model), which is the Abelian HQET. In the Abelian case, this non- 
renormalization holds to all orders due to the exponentiation theorem. In HQET. this is only 
true at one loop. 

Now we shall discuss the renormalization of g. Due to the gauge invariance, all 5's in 
the Lagrangian are equal. The coupling of the HQET quark field to giuon is thus identical 
to the usual QCD coupling, where the HQET heavy flavour is not counted in the number of 
flavours Uf. In order to find Za we need to renormalize the heavy-quark - gluon vertex and all 
propagators attached to it. We have already calculated Zq. The renormalization of the giuon 
propagator is well known (Fig. [6]): 



Za = 1 



Ca 
2 



13 

y 



-TFUf 



as 



(2.32) 
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Figure 6: One-loop gluon self-energy 



Let's introduce the vertex 




(2.33) 



where starts from one loop. When expressed via renormalized quantities, the vertex should 
be r = ZrTr, where the renormalized vertex is finite at e — > 0. A physical matrix element is 
obtained from the corresponding vertex by multiplying it by the wave-function rcnormalization 

1/2 

constant for each external leg. In our case, 

gotZqZ]/^ = .gf rZy^ZrZgZy^ = finite. 



Therefore, Zl/^ ZyZqZ^'' must be finite. But the only minimal (|2.24p rcnormalization constant 
finite at e ^ is 1: 

Z^ = {ZrZQy^Z^K (2.34) 

At one loop the HQET vertex is given by two diagrams (Fig. [7]). It is very easy to calculate 
the first one. It contains two heavy denominators which can be replaced by a diff'erence: 



A/2 



1 



1 



1 



1 



(ko + Lu)[kQ + llj') lu' — lu \kQ + uj ko + uj' 
We get the difference of the self-energies: 



A^ 



^ _ Ca \ SKl^EH . 
2CfJ uj' -uj 



(2.35) 



(2.36) 



This result can also be obtained from the Ward identity. The UV divergence of this contribution 
is 

Ca^ 



A'/ 



(a -3) 



(2.37) 




Figure 7: One-loop HQET vertex 
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The second diagram is more difficult. It has been calculated in [5]. Now we only need 



its UV divergence, and it should be 



We may nullify all external momenta. After 



that, the diagram will contain no scale and hence vanish. It will contain both UV and IR 
divergences which cancel. Therefore, we'll have to introduce some IR regularization to get the 
UV divergence. We have 



Ca 



f^gov 



k ■ V 



fc2 



T^i9PP'-c'^)mV-^'^{k,-k,0), 



(2.38) 



where the three-gluon vertex is 

y^i^^'^Mpi,P2,P3) = (P3 -P2)^^5^^^^ + {pi-PsT^g^'^' + {P2-Pir'9^'^^ ■ (2.39) 

It vanishes when contracted with the same vector in all three indices: 



0, 



and when contracted in two indices with the corresponding momenta: 

y/^lM2M3(p^^p2,P3)PlMlP2^. =0. 

Therefore, ^'^ and terms vanish: 

d''k k^ -{k- w)2 



Averaging over k directions (fc • u)^ — > k^/d, we get 



A^«,, = iCAgU [ 1 - 



d'^k 1 

(2^(fc2)^ 



The UV divergence of this integral can be obtained by introducing any IR regularization, e.g., 
an IR cut-off in the Euclidean momentum integral or a small mass: 



/ 



d'^k 1 



{27T)d (fc2)2 



UV 



We arrive at the UV divergence of the second vertex diagram: 



:CA{l~a) 



47re 



From (p:^ and (piiT]) we obtain 



Zr 



CFia-3) + CA- 



Aire 



The product which appears in (j2.34p is 



ZrZq = 1 + Ca 



a + 3 tts 
4 Ane 



(2.40) 



(2.41) 



(2.42) 



(2.43) 
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In the Abelian case ZtZq = 1 to all orders, due to the Ward identity. This is why we only got 
the non-abelian colour structure Ca in (|2.43p . Finally, combining this with (|2.32p . we see that 
the a dependence cancels, and 

^a = l-/3o^, fio^^4^A-\TFnf. (2.44) 
47re 3 3 

Thus we have derived the one-loop (3 function of QCD {nf does not include the HQET heavy 
flavour Q). Most textbooks use the massless-quark - gluon vertex or the ghost - gluon one (see, 
G-g-i l2])- In the later case, calculations are a little shorter. The HQET derivation presented 
here is as short as the ghost one. 



3 a parametrization 
3.1 General formulae 

a parametrization of Fcynman integrals (including those containing numerators) is discussed 
in many textbooks, see e.g. [5]. Here we shall discuss HQET integrals; all rules can be trivially 
obtained from [S], though they were not yet explicitly stated in the literature. 

First let's calculate the one-loop diagram (|2.9p (Fig. [T|) using a parametrization 



1 



daa'^-'e-'"^ . (3.1) 



a" r(n) Jo 
We get 

^, , I da dp 13'"'-^ d'^k , ak^ + 2/3(fc + p) 

T{ni)T[n2) J 

We shift the integration momentum 

k = k' — —V 
a 

to eliminate the linear term in the exponent: 



fl2 

X = ak'"^ ~ — + 2(3uj. 
a 

Now it is easy to calculate the momentum integral: 

d'^k e"''' = i f d'^ks e""''^^ = i (^) '^^^ . (3.2) 



Therefore, 

(-2^)'^-"^-'"^/(7ii,n2) ^ ^, l^, , / daa"=-id/3/3"i-ia"'^/2expf-^ + 2/3cj^ 

r(?ii)r(n2) J V " y 

Now we make the substitution [3 = ay and integrate in a: 

1 d/2 — ni — n2 



\, ' / dyy-^''[y{y-2u;)y 

r(7ii)r(n2) Jo 
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The HQET Feynman parameter y has the dimensionaHty of energy and varies from to oo. 
The y integral can be easily calculated in T functions, and we again obtain (|2.17p . 

Now we shall consider the most general HQET Feynman integral without numerators. Any 
HQET diagram contains a single heavy line and has the form 

nd^ki 1 

where ki are loop momenta G 

La = ml - ql- iO , He = -2(jc ■ v - iO 
are light denominators (a, b € [1, Ni]) and heavy ones (c, d € [1, Nh]): 

ga = ^ Naiki + ^ NanPn , 9c = ^ Ndki + ^ N^nPn , 

Pn are external momenta {n,m € [1, — 1]), and the coefficients express momenta of 
propagators via the loop and external momenta (these coefficients are equal to or ±1). Using 
the a representation (j3.ip for all lines, we obtain 



I 



d'^h X 



X = Y,aa{ql-ml) + 2Y,Pcqc-v. (3.4) 
Dimensionalities of the parameters are aa ~ lA^^, ^ 1/M. The exponent is 

X^Yl • fc, - 2 ^ • h + Y , (3.5) 

where 

My = aaNa^NaJ , 

Qi = - ^ CtaNaiNanPn ~ PcNd , 

Now we shift the loop momenta ki = k[ + Ki to eliminate linear terms: 
Then 

X = Y, K ■k'.-Yl ^^'U 'Q^■Q3+y■ (3-8) 
Performing the Wick rotation to Euclidean k^ and integration in the loop momenta, we obtain 

A{a) + Ai{a, P) + A2{a, f3) 



X exp 



(3.9) 
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where 



D{a) ^ dot A/, 
D{a) 
D{a) 



= ^ ^aaabNaiNanNbjNbmPn ' Pm - ^ OlaNanNamPn ' Pm , 

2 M~J^aaf3cNa^NanNcjPn ' " 2 ^ /ScNcnPn ' ?^ , 



= E Vc/3diVc,^d, . (3.10) 

The polynomials D{a), A{a), Ai{a,P), A2{a,f3) have dimensionality 1/M^^. The fmiction 
D{a) is homogeneous in aa of degree L. The function ^(a) is homogeneous in aa of degree 
L + 1 and linear in p„ • pm- The function ^i(a, (3) is linear in /3c, of degree L in a^, and linear 
in Pn ■ V. The function ^2(0?, /?) is quadratic in (3^ and of degree L — 1 in a^; it docs not contain 
momenta. 

It is always possible to calculate (at least) one integration in p.9[) . Let's insert S (^ aa — i]) drj 
under the integral sign, and make the substitution aa = ?ya;aj Pc ~ VUc- Then the rj integral is 
a r function: 



r(Ena + E^c-xf) f 

mMunnc) J 



UdXaX^^-^ Udycy^'-'S{j:Xa-l) 



D{x) 

(3.11) 

The ordinary Feynman parameters Xa are dimensionless and vary from to 1; the HQET 
Feynman parameters j/c have dimensionality of energy and vary from to 00. 



3.2 Graph-theoretical rules 

The polynomials D{a), ^(a), Ai{a,f3), A2{a,P) can be extracted directly from the diagram, 
see [6]. Wc shall formulate the rules and illustrate them by an example shown in Fig.[8l 







' P2 






H 




/ 




Pi 


1 


2 



Figure 8: An HQET vertex diagram 

1. Cut a few light lines so as to get a connected tree, form the product of aa of the cut 
lines. D(a) is the sum of all such products. 




D ~ a\a2 + a\a-i + 0203 . 
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2. Cut a few light lines to get two connected trees (the heavy line is in one of the two parts). 
Form the product of aa of the cut lines and multiply it by (— P^), where P is the momentum 
flowing from one connected part to the other one. A{a) is the sum of all such terms. 



3. Cut a single heavy line and a few light ones to get two connected trees (now the heavy 
line enters one connected part and leaves the other one). Form the product of /3c of the cut 
heavy line and aa of the cut light ones and multiply by (— 2P • w), where P is the momentum 
flowing from the first connected part to the second one. Ai(a,/3) is the sum of all such terms. 






A\ = -2pi ■ V (3iai{a2 + aj,) - 2{pi + P2) ■ v (3ia2a3 
- 2{pi +P2) ■ V (32a2{ai + a^) - 2pi ■ v f32aia3 . 

4. Cut a few light lines to get a connected diagram with a single loop in such a way that 
this loop contains at least one heavy line. Sum Pc of the heavy lines belonging to the loop, 
square the sum, and multiply by aa of the cut light lines. Sum all terms. 






A2 



a2f3f + aaiPi + P2? 



These rules can be simplified a little. Suppose pi is the residual momentum of the incoming 
heavy line, and let's route it along the heavy line. Then the exponent X (|3.4[) contains 2pi ■ 
V (3c. We can add this expression to the exponent in (j3.9p . and then set pi = while 
calculating Ai {a, (3). 

There is an analogy between Feynman diagrams in a representation and electrical circuits 
(Table [1]). The average momentum flowing through a propagator corresponds to current. The 
flrst Kirchhoff rule is satisfled: the sum of momenta flowing into a vertex vanishes. Light lines 
are resistors aa, and heavy lines — voltage sources (3cV (batteries with zero internal resistance, 
the voltage does not depend on the current). The second Kirchhoff rule says that the sum of 
voltages along a loop (say, loop i) must vanish. These equations are nothing but the equations 
^MijKj = Qi which determine the average loop momenta Ki (see p.7p ). 



The Joule heat ^ aa'Za pl^^^ 2 times the energy consumption by the voltage sources ^ (3cqc 



gives 



A{a) + Ai{a,l3) + A2{a,[3) 
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Current 


Voltage 


• > • 

a 






c 







Table 1: Analogy with electrical circuits 



The case ao ^ or /3c ^ corresponds to a short circuit (line shrinks to a point); the case 
Qfa — > c>o — no contact (the line is removed). 

Generalization to integrals with numerators is straightforward |6j . Suppose we have a poly- 
nomial V{qa, Qc) inserted into the numerator of (|3.3p . Then we can add 

2 ^ go • + 2 ^ • f?c 
to the exponent X (|3.4p . and apply the differential operator 

'1 d 1 d 



V 



2 ' 2 dljc 



to the result at = 0, 77c = 0. Before this step, the only difference is the substitution 
PcV — > PcV + 77c, and the fact that a light line a can be also considered "heavy" in Ai and A2 
calculations, with playing the role of (iv. Let's formulate the rules to calculate Ai, A2. 

3'. Cut a single heavy line (say, c) and a few light ones to get two connected trees, and form 
the product of —2{(3cV + i]c) ■ P, where P is the momentum flowing from the first connected 
part to the second one; multiply this product by aa of all cut light lines. Or cut a single light 
line (say, a) and a few light ones to get two connected trees, form the product — 2^^ ■ P, and 
multiply it by at of these additional cut lines. Here the first and the second connected parts 
are defined by the direction of the momentum of the first cut line {qc or qa for a heavy or light 
line). Ai(a, (3, ^, rj) is the sum of all such terms. 

4'. Cut a few light lines to get a connected diagram with a single loop. Sum (3cV + r]c or 
of the heavy or light lines belonging to the loop, square the sum, and multiply by aa of the cut 
hght lines. Sum all terms to get v42(a, /3, ^, 77). 



4 HQET propagator diagrams 

4.1 Two loops 
4.1.1 Diagram 1 

We have calculated the one- loop HQET propagator diagram by three different methods (Sects. [272 
and lS.ip . Now we shall consider two-loop propagator diagrams in HQET. There are two generic 
topologies of such diagrams (Fig. [H]). This means that all other possible topologies can be ob- 
tained from these ones by shrinking some lines. The method of calculation of these diagrams 
has been constructed in [7]. 

The first diagram (Fig. [T0|) is 

- ^ J orm'^^-Bp ^ (-2-)'^-— . (") 
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Figure 9: Generic topologies of two- loop propagator diagrams 




Figure 10: Diagram 1 



where 

Di ^ -2{kw + Lu) , D2 = -2{k2o + uj) , 
D3 = ~kl, Di = -kl, D^^-ik^-k^f 

(the power of —2lj is fixed by dimensionality). It is symmetric with respect to (1 <-> 2, 3 <-> 4), 
and vanishes if two adjacent indices are < 0. 

If 715 = 0, the diagram is the product of two one-loop ones: 

713 "4 



/(ni,n2,n3,«4,0) 




/(?ii, 713)7(712,714). (4.2) 



71l 7l2 

If 7ii = 0, wc first calculate the inner masslcss loop. The one-loop massless diagram is 
1 [ d'^k 



_ ('_„2Nd/2-ni-r 

ind/^ J [-{k+pf -iOf'[-k^ -iO]""' ^ ^' 

^, . r(-d/2 + 7ii + 7i2)r(rf/2-7ii)r(d/2-7i2) 

G(71l,?7.2) = 



•G(71l,7l2) : 



r(7ii)r(7i2)r(d - Til - 712) 

This gives the coefficient G{n3, 715), and shifts the power 714 by 713 + n<^ ~ d/2: 



(4.3) 




"3 



^^/^^^'^^ 714 713 -I- 715 - d/2 




n2 



715 



1(0,712, 713,714, 715) = G{n3,n5)I{n2,n4 + 713 + 715 - d/2) 



(4.4) 
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(the case 712 = is symmetric). If na = 0, we first calculate the inner HQET loop. This gives 
the coefficient /(ni, 715), and shifts the power n2 by ni + 2n5 — d: 




7I5 713 





77l 772 + ni + 2775 ^ d 

Til n-2 

/(711, 772, 0, ^^4, f^5) = I{ni,n5)I{n2 + 77i + 2775 - d, 714) (4.5) 

(the case 7t4 = is symmetric). 

But what can wc do if all 5 powers of denominators are positive? We shall use integration 
by parts [S]. Integral of any full derivative over the whole space of loop momenta is zero. When 
applied to the integrand of (|4.ip . the derivative 

A^^2. + ^2fc2 + ^2(fc.-M. 
dfc2 D2 Di D5 

Applying {d/dk2) ■ k2 or {d/dk2) ■ (fc2 — fci) to the integrand, we obtain zero integral. On 
the other hand, we can calculate these derivatives explicitly. Using 2fc2 ■ v = —D2 — 2uj, 
2(fc2 — fci) • fc2 = 1?3 — D4 — -D5, we see that applying these differential operators is equivalent 
to inserting 

d - 7i2 - 775 - 2774 - 2w-^ + ^{Ds - Di) , 

d - 772 - 774 - 2775 + + ^{Da - D5) 

L>2 Ui 

under the integral sign. These combinations of integrals vanish. These recurrence relations are 
usually written as 

[d- 772 -775 - 2774 + 7722+ +7755+ (3" -4")] 7 = 0, (4.6) 

[d - 772 - 774 - 2775 + 7722+1" + 7744+(3" - 5")] / = , (4.7) 

where, for example, 1~ lowers 774 by 1 and 2+ raises 772 by 1. Applying {d/dk2) ■ v, we obtain 
a (less useful) relation 

[-27722+ + 7744+(2- - 1) + 7755+ (2- - 1")] 1 = 0. (4.8) 

A useful relation can be obtained from homogeneity of the integral (j4.ip in uj. Applying 
Lo^d/dio), we get the same integral times its dimensionality 2(d — 773 — 774 — 775) — rii — ?72. On 
the other hand, we can calculate the derivative explicitly: 

[2(d - 773 - 774 - 7I5) - 771 - 772 + 771 1+ + "-22+] 7 = 0. (4.9) 



This homogeneity relation is not independent: it is the sum of the (d/dk2) ■ fc2 relation (|4.6p 
and its mirror-symmetric (d/dki) ■ ki one. 

A particularly useful relation can be obtained by subtracting the 1~ shifted homogeneity 
relation from the {d/dk2) ■ {k2 - fci) relation (|t?7)) : 

[d — 77l — 772 — 774 — 2775 + 1 ^ (2(d — 773 — — 775) — 771 — 772 + l) 1^ 
+ 7744+(3" - 5")]/ = 0. (4.10) 
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Solving it for the / with the unshifted indices, we obtain an expression for /(ni, n2, n3, n4, ns) 
via 3 integrals: 



(2((i - 713 - n4 - ns) - m - 712 + 1)1 +7i44+(5 )j 



ni 



712 - 714 — 2715 + 1 



(4.11) 



Each of them has tii + 713 + 715 reduced by 1. Each application of (|4.1ip moves us closer to the 
origin (Fig. Ilip . Therefore, after a finite number of steps, any integral /(tii, 772, 773, 774, 715) will 
be reduced to the trivial cases in which one of the indices vanishes. 



Figure 11: A single step of the integration- by-parts reduction 



4.1.2 Diagram 2 

The second diagram (Fig. [T2|l is 

1 f d'^ki d'% 



where 



= (-2cc>)2'^-"l-"=-"-^-2("^+"=) J(77i, 772, 773, 774, ^5) , (4.12) 



Di = -2{kw + uj) , D2 ^ -2{k2o + Lj) , D3 = -2(fcio + fc20+t^), 
D4 = -kj , D5 ^ -kl 

(the power of — 2cj is fixed by dimensionality). It is symmetric with respect to (1 <-*■ 2,4 ^ 5), 
and vanishes if 774 < or 775 < or two adjacent rii 3 are < 0. 

This integral is trivial if 773 = or 771 2 = 0. In general, it has 3 linear denominators and 
only 2 loop momenta; therefore, these denominators are linearly dependent: 

Di + D2 - D:i = -2uj . (4.13) 

Inserting this combination under the integral sign, we obtain 

J= (1" + 2- -3")J. (4.14) 
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Figure 12: Diagram 2 



Each application of this recurrence relation reduces nl+r^2 + ?^3 by 1. Therefore, after a number 
of such steps any integral will reduce to the trivial cases (Fig. [TT|) . 

The integral (|4.12p can contain a power of ki ■ k2 in the numerator; this scalar product 
cannot be expressed via the denominators. However, this is not a serious problem [9j. 

Let's summarize. All scalar integrals belonging to the two generic topologies of Fig. [51 with 
any indices rii (and with any power of fci • /c2 in the numerator of the second integral) can be 
reduced to linear combinations of two master integrals 



with coefficients being rational functions of d. Here the n-loop HQET sunset integral is 



This reduction can be done using integration by parts [7] (see also [5]). 

4.2 Three loops 
4.2.1 Reduction 

There are 10 generic topologies of three- loop HQET propagator diagrams (Fig. [T^. 

All these integrals, with any powers of denominators and irreducible numerators, can be 




(4.15) 




(4.16) 
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Figure 13: Generic topologies of three-loop propagator diagrams 
reduced [S] to 8 master integrals: 

^^^^^^/3, (4.17) 

^^'''^'^^^^^ - ^1,1,, (4.18) 





= h , (4.19) 



/^"""h/^^ _ 1^1(6 - 2d, 1) , _ 3d - 7 /J 



/2/(5-2d, 1) 2d-bh 



^3 = 77^^ ^^3, (4.20) 



/ff'^V|\ _ G'f/(l,4-d)^ _ ,3d-7G^ 



^3 = -2-^^ 7^/3, (4.21) 



G2/(l,3-d) " d-3 G2 





(4.24) 



using integration by parts. Here the n-loop HQET sunset /„ is defined by (|4.16p . and the n-loop 
massless sunset is 



Gn 



r(l + n£)r"+i(l -e) 



{n + 1- n|)„ ((n + l)f - 2n ~ l)^ r(l - (n + 1)£) 



(4.25) 

This reduction algorithm has been implemented as a REDUCE package Grinder [5]. It is analo- 
gous to the massless package Mincer [10]. The first 5 master integrals can be easily expressed 
via r functions, exactly in d dimensions. The next two ones reduce to two-loop ones with a 
single e-depcndent indes0 (Sects. IT2.2l and B:.2.3p . The last one is truly three-loop (Sect. ITTl)) . 



4.2.2 J(l, 1,71,1,1) 



Here we shall calculate the integral J (|4.12p (Fig. for arbitrary powers of denominators. To 
this end, we shall first consider the one-loop diagram with two different residual energies toi 
and LU2 (Fig. [Hi): 

where 



Di = -2(fco + c^i), Da = -2(fco + c^2), Dj, = ^k^ . 

If ni^2 are integer, this integral can be easily calculated by partial fraction decomposition 
(Sect.[2J|). 




UJl fco +071 fco + UJ2 0J2 




Figure 14: One-loop diagram 
Closely following Sect. 12.21 we first integrate in d'^^^k: 

d-l\ ^ (1,2 ^(<i-l)/2->^3, 



(fc 



"eg I 



'dk 



^i/2r(n3) J {-2uji - 2ikEor" {-2uJ2 - 2ikEoy 

. r(n3-^) 

! — = — - CO 

7ri/2r(n3) 



■^dfc 



(2fc - 2wi)"i(2fc - 2w2)"2 



^Grinder uses B4 = l3l'f/I2 and S5 = /3GJ/G2 as elements of its basis instead of 114. 20t and 1 14. 211 . 
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We obtain [TT] 



/ = I{ni + n2, 713)2^1 



rii,ni + ^2 + 2^3 - d 
ni + 712 



1 - — ) (-2^2)''""'""'"^"' • (4.27) 

LJ2 



One can easily check that this resuh is symmetric with respect to {lji ^ 1-^2, ni <-> 0.2), using 
properties of hypergeometric function. If uji ~ 102-, it reduces to l(n\ + n2, 713). 
Let's also calculate this integral using a representation (Sect. [3TT|) : 



J 



r(7ii)r(772)r(7i3) 



(/3i + /32)' 



X exp 



+ 2(wi/3i +W2/32) 



Now we make the substitution /3i^2 = ctyi,2 and integrate in a: 

"2/1 2/1 "2/2 2/2' [(2/1+2/2) - 2(cjiyi + t^22/2)J 
After the substitution yi = yx, 7/2 ~ 2/(1 ^ 2^), the integral in y can be taken: 



^_ r(r7l +712 +773 - f) 

r(7ii)r(7i2)r(7i3) 



^ _ r(f - 7i3)r(7ii + 712 + 2713 - f ) 
r(7ii)r(7i2)r(773) 



dxx'''-^ (1 - a;)"^-i [-2t^ia; - 2u;2(l - x)Y 



-rii— n2— 2n3 



(4.28) 



And we again obtain (j427l) . 

Finally, we shall derive the same result in coordinate space fFig. [Mb. see Sect. 12.2] ): 



r(i-7i3) 



4r(7ii)r(7i2)r(7i3) 



dtidt2 e*('^i*i+'^2*2) 



ni —1 



1^2 

T 



71-2—1 



7(^1+^2) 



2Tl3-dl 



The substitution = tx, t2 = t{l — x) reduces this expression to (|4.28p . 

Now we return to our main problem — calculating J = J{ni, 712, 713, 714, 715) (Fig. I12p with 
arbitrary indices. We set — 2w = 1; the power of — 2w can be reconstructed by dimensionality. 
Substituting the one- loop subdiagram (|4.27p . we have 



J 



/(711 +713,714) 

j_-d/2 



d<^k 



(_fc2)n5(l _ 2fco)"^ 

^ Q _ 2A: )''-"i-"3-2n4 p f 711,711 + 713 + 2714 - d 
° ^ ^ V 711+773 



-2fco 



1 - 2k^ 



/(ill +7l3,7l4)r(n5 



+00 



(fc|o)('^-l)/2-"=dfc 



BO 



X 2F1 



7r''/2r(7i5) 

71l, Til + 7I3 + 27I4 — d 
111 + 7I3 



(1 - 2iA:£;o)"i+""+"''+^"-'~'' 

-2lfc£;o 



1 - 2ik 



EO 
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We can deform the integration contour (Fig.[3l ksQ = iz/2): 



/(ni +n3,n4)r(n5 - ^) 
J ~ ^r-T. ^ — rh. — -, r^-^ COS 



2d-2n5-%'i/2r(n5) 

^d_2„5-l^^ / ^^^^^ _|_ _,_ - d 

711+713 



(^2 + I)ni+n2+n3+2n4-d ^ 

Now we substitute the series 



z + 1 



2F\ 



a, b 
c 



_ r{c)r{b) ^ T{n + a)T{n + b) „ 

r(a) A.r(77 + i)r(77 + c)'' 



and integrate term by term. The resuh is 

I{ni + 713, ni)T{n5 - '^)T{d - 2n^)r{ni +712 + 713 + 2714 + 2715 - 2d) 



J 



2d-2«5-l7rrf/2r(n5)r(71i + 712 + 713 + 2714 - d) 

Til , 7ii + 713 + 2714 — d,d ~ 2715 _^ 



COS 



TT 7I5 



71l + 7I3, 7ll + 7l2 + 7I3 + 2?l4 - d 

Using p.l6p . we can simphfy this resuh: 
7(711,712,713,714,715) = 

r(| - "4)r(f - 7i5)r(7ii + 713 + 2714 - d)r(7ii + 712 + 713 + 2714 + 2715 - 2d) 



X 3F2 



r(7i4)r(n5)r(7ii + n3)r(7ii + 712 + 713 + 2714 — d) 

7ii , 7ii + 713 + 2714 — d,d — 2715 

7ll + 713, 71l + 7l2 + 7I3 + 2714 - d 



(4.29) 



It was first derived in coordinate space [9l [3] . Checking the symmetry (tii <-> 7i2 , 714 
requires using some 3i^2 identities. 



4.2.3 /(1, 1,1, 1,71) 

This diagram has been calculated in [12j using Gegenbauer polynomial technique in coordinate 
space [13]: 




7(1,1,1,1,71) 



2r(27i-7 + 3)r(27i- 2d + 6) ^ 
(7i-d + 3)r(37i-2d + 6) ^ ^ 

-r(d-7i-2)r2(7i-(i + 3) 



r(i-i)r(|-7i-i) 
r(fi- 2) 

71 - + 3, 71 - 7 + 3, 271 - 27 + 6 
71 - 7 + 4, 371 - 27 + 6 



(4.30) 



Some details of this method are discussed in [3]. 
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4.2.4 Inversion 

The last three-loop master integral has been calculated [TS] using inversion. We shall first 
consider inversion relations at one and two loops. The one- loop massive on-shell integral defined 

by 



[m2 -{k + TOU)2 - iOf' [-k^ - iO] 

can be written in terms of the dimensionless Euclidean momentum K — /m 



I 



{K^ - 2iA'o)"i(iir2)n2 

Similarly, the one-loop HQET propagator integral ()2.9p expressed via K = kE/{—2uj) is 

^ ^ , — -- — = Tr'^/^/fni, 712) . 

(1 - 2ii^o)"i(A'2)"2 ^ 1' ^' 

Inversion K = K'/K'^ transforms the massive on-shell denominator into the HQET one: 

2 o... _ l-2*i^^ 



- 2iKo 



K 



n2 



Therefore, 

n2 d — rii — 712 





ni ni 



M( \ T( ^ A r(rf-ni-2n2)r(-rf/2 + ni+n2) , . 

M[ni,n2) = I(ni,d- ni - n2) ^ — — ——: r . (4.32) 

r(7ii)r(d - rii - 1X2) 

Similarly, at two loops we obtain |14j 









d 


















ni 


n2 









\\d- 




- nspT < 






ni 


n2 





(4.33) 



This relation is less useful, because the HQET diagram in the right-hand side contains two 
non-integer indices. 

At three loops we have [3j 

d — ne — nj — ng 
d — ni — n^Sn-. — 712 — 714 

"-5 — n-jL^ng, \— - Tig 

>n5 




ni 




"2 




d-n^- 




nj - 




ni — n4 Jl 






- "-2 - «5 


- rtg j-T < 


>7i6 nj' 




^ "-7 



(4.34) 



?Zl 713 ?12 



ni 713 ?^2 
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In particular, the HQET ladder diagram with all indices = 1 is convergent; its value at d = 4 
is related [15] to a massive on-shell diagram 



-5C5 + I2C2C3 • (4.35) 



by the second inversion relation. This is one of the on-shell three-loop master integrals, and its 
value at d = 4 is known [16|, I17j. Calculating this ladder diagram with Grinder: 



(d-3)^ 3 136 (d-3)(2d-5)(2d-7) ^ 



2 (3d - 7)(3d - 8)(81d3 - 891^^ + 3266d - 3988) 



(d-4)2 ^ 3 (d-4) 



2 



(d-4)4(2d-7) 

jd- 3)(3d - 7)(3d - 8)(3d - 10)(3d - 11) If (d - 3)(3d - 7)(3d - 11) G? 
^ (d-4)3(2d-5)(2d-7) /2 (d-4)3 G2 ^ 

and solving for the most difficult HQET three- loop master integral Bs (|4.24p . we obtain the 
£ expansion of this integral up to 0{e). This concludes the investigation of three- loop master 
integrals, and allows one to solve three-loop propagator problems in HQET up to terms 0(1). 



4.2.5 Applications 

Using this technique, the HQET heavy-quark propagator has been calculated up to three 
loops jl8J, and the heavy-quark field anomalous dimension (obtained earlier by a completely 
different method [17j ) has been confirmed. The anomalous dimension of the HQET heavy-light 
quark current has been calculated [18j . The correlator of two heavy-light currents has been 
found, up to three loops, including light-quark mass corrections of order m and |19| . The 
quark-condensate contribution to this correlator has been also calculated up to three loops [19j . 
Its ultraviolet divergence yields the difference of twice the anomalous dimension of the heavy- 
quark current and the that of the quark condensate, thus providing a completely independent 
confirmation of the result obtained in |18| . The gluon-condensate contribution has been calcu- 
lated up to two loops [T9j (at one loop it vanishes). 



5 On-shell HQET propagator diagrams with mass 
5.1 Two loops 

On-shell HQET propagator diagrams vanish if all flavours (except the HQET one) are considered 
massless, because loop integrals contain no scale. If there is a massive flavour (c in the 6-quark 
HQET), such diagrams are non-zero. They first appear at two loops. They are used, e.g., to 
calculate on-shell renormalization constants in HQET. 

Let's first consider [2^ a class of such integrals (Fig. [T5|) 

/f(k^)d'^k 
j^rjjjn, , where D, ^ -2k ■ v - lO , D2 ^ -k^ - iO , (5.1) 
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Figure 15: Diagram with a single HQET line 



and f{k'^) is an arbitrary function. We can construct an identity in which f'{k'^) terms cancel: 



^■(k-2£lv]J^ 

dk y ^r^^r 



fi-m -2-4(711 + 1):^ 
^1 



(5.2) 



Integrating it, we obtain an integration- by-parts relation 



(d - ni - 2)F(ni, 712) = Mm + l)l++2- F{ni,n2) . 



(5.3) 



Let's call integrals with even ni apparently even, and with odd rii — apparently odd (they 
would be even and odd in v if we neglected iO in the denominator) . These two classes of integrals 
are not mixed by the recurrence relation (j5.3p . We can use this relation to reduce all apparently 
even integrals to vacuum integrals with ni ~ (Fig. [16]). Apparently odd integrals with ni < 
can be reduced to ni = — 1. Substituting = — 1 to (|5.3p . we see that these integrals vanish, 
and hence all integrals with odd ni < vanish too. Apparently odd integrals with ni > can 
be reduced to ni = 1 (Fig. [TO]) : however, they are not related to those with ni = —1. 



Ari2 




Figure 16: Recurrence relation 
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The solution of the recurrence relation can thus be written as 

-F 0,722 + - 



(-4)- 
2i-»i 




-ni/2_ (2) 



r(^) r(i) 
r 



r(^)r(^) 



ni — 1 



even ni, 

odd Til > 0, 
odd ni < 0. 



(5.4) 



Some of these properties can be understood more directly. If ni < 0, iO in Df"^ can be 
safely neglected; averaging this factor over k directions, we obtain for odd ni and the upper 
formula in (|5.4p for even ni. It was suggested [21] that this last formula can also be used for 
even ni > 0, but the proof (presented here) only appeared in [20j . 




Figure 17: Two-loop diagram 
Now let's consider the two-loop diagram fFig. Wf} 

where 

£>! = -2fci • w - iO, D2 = -kf-i0, 

D3 = l-kl-iQ, D4 = 1 - (fci -I- fca)^ - iO. 

It is symmetric with respect to 3 <-> 4, and vanishes if 713 or 714 is integer and non-positive. It 
can be calculated using a parametrization |20j : 

F(71l,7l2,7i3,7l4) = (5.6) 

- "2)r(^ + nz + n3)T{^ +n2 + ni)r{^ + 712 + n3 + 7t4 - rf) 

2r(7ii)r(n3)r(?24)r(^)r(7ii + 2^2 + n3 + n4-d) 

In full accordance with (j5.4|) . integrals F{ni,n2,n3,n4) with even ni reduce to F{0,n2 + 
711/2,713,714) (this is a well-known two-loop vacuum integral [22]); those with odd ?ii > 
reduce to F{1, 712 -I- {ni — l)/2, 773, 714); and those with odd Tii < vanish. All apparently even 
integrals are proportional to the single master integral 
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and apparently odd ones — to 



Jn = 




^4d-9 2 r(5-2rf) 



r2(2-d/2) 



(5.8) 



Integrals (|5.5[) can also contain powers of {2k2 + ki) ■ v in the numerator; see [20j for details of 
their evaluation. 



5.2 Three loops 
5.2.1 Reduction 

There are two generic topologies of three-loop on-shell HQET propagator diagrams with a 
massive loop fFig. 118(1. Algorithms of their reduction to master integrals, using integration by 
parts identities, have been constructed [50] by Grobner bases technique 





Figure 18: Topologies of three-loop on-shcU HQET propagator diagrams with mass 
All apparently even integrals of the first topology reduce to 







while apparently odd ones to 



^ ^ 



All apparently even integrals of the second topology reduce to 




while apparently odd ones to 
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The master integrals 




can be easily expressed via F functions. The master integral 



has been investigated in detail 



5.2.2 A master integral 

Now we shall discuss the integrals 



n2 n2 




111 ni 



(5.9) 



(/ill is one of the master integrals). Several approaches have been tried [20l[26]. The best 
result was obtained using a method similar to [S^ . 

First we consider (following [55]) the one-loop subdiagram 




1 



dko (F^'^k 



ni iTT'^/z J [-2{ko + po) - «0]"i - iO]"2 

After the Wick rotation, we integrate in d'^^^k: 

T{n2- id -l)/2) (fc|o + l)('^-^)/^-"- 

'"^"^ - ^V2r(n,) i_ '^'^^ i-2po - 2^fc,o)- ■ 

If Po < 0; we can deform the integration contour (Fig.[T 



, , F(n2 - (rf- l)/2) 



7ri/2r(,i2) 



TT ?7.2 

2 



dk 



(fc2 _ iYd-l)/2-n2 

{2k - 2po)"i 



This integral is 



Inin2 {Po) 



r{ni + na - 2 + £)F(ni + 2^2 - 4 + 2e) 

r(n2)F(2(ni +n2 -2 + e)) 

/ ni,ni + 2n2 -4 + 2e 1 , , 
^2Fi[ 3 :^(l+Po) 



(5.10) 
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or, after using a 2-F1 identity, 



r(ni + na - 2 + £)r(ni + 2n2 - 4 + 2e) 
r(n2)r(2(ni+n2-2 + e)) 

^ TP f 5"!' + '^2 - 2 + £ 

X 2^^! „ I „ 3 , . 

??1 + ^2 — 2 + ^ 



This result was obtained |26j using the HQET Feynman parametrization: 

r(ni+n2-2 + £) r n.-l^,2 0..„^n2-n,-„.-e 



Inin2{Po) 



dy. 



r(ni)r(n2) Jo 

This integral at po < gives (|5.11l) (a similar expression has been derived in [27]) 



(5.11) 



C 



Wo 



Figure 19: Integration contour 

Now we can integrate in d'^~^p in the three-loop diagram: 
r(n3-3/2 + e) '•+°° 



^i/2r(n3) 



2 \3/2-n3-e 



dpEO ■ 



(5.12) 



The square of 2^1 in (|5.1ip can be expressed via an 3F2 using the Clausen identity. We 
analytically continue this 3F2 from 1 + p%q > 1 to 2; = 1/(1 + p'eq) < 1 and integrate (|5.12p 
term by term. The result contains, in general, three 4^3 of unit argument. 
A convergent integral /122 is related to the master integral /m by 



(d-3)^(d-4)(3d-8)(3d-10) 
8(3(i- ll)(3d- 13) 



hi 



For this integral, we obtain 

-fl22 _ 1 

" "2^ 



r3(i + £) 



"4^3 



1, i - e, 1 + £, -2£ 

I + £, 1 - £, 1 - 2£ 

r2(i-£)r3(i + 2£: 



1 + 4£ r2(i + £)r(i - 2£)r(i + 4£) 

1 r2(i-£)r''(i + 2£)r(i-2£)r2(i + 3£) 
1 + 6£ r4(i + £)r(i + 4£)r(i - 4£)r(i + 6£) 



i,l + 2£, -£ 
I +2£,1 -£ 



(5.13) 
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Expansion of this result up to agrees with [26] 

/i22 7r2r3(l + 2e)r2(l + 3£ 



r3(i + e) 3 r6(i + £)r(2 + 6£) 



(5.14) 



This equality has also been checked by high precision numerical calculations at some finite e 
values. This conjectured hypergeometric identity can also be rewritten in a nice form |28| 



232(£)3i^2 



i,l + 2£,-£ 



2£,1 



1 +53(£) = 0, (5.15) 



where 



6"(£ 



r(i + £) 



6(ri£)(l + 2n£) 



We have no analytical proof. 



5.2.3 Other master integrals 

Other master integrals were calculated |2D] using Mellin-Barncs representation (sec, e.g., U). 
Now we shall discuss a simple example of this technique. Let's consider the one-loop propagator 
diagram with two massive lines: 



1 



Using Fcynman parametrization, 
V{n\ + 712) 1 



r(ni)r(n2) i7r''/2 
V(nx+n2) 1 



r(ni)r(n2) i7r'^/2 
After the shift k — k' — xp: 

r(ni+n2) 1 



[m2 - fc2]«i [to2 _ (fc +p)2]n2 ' 

da;x"=-i(l -a;)"i-id'^fc 
[(to2 - A:2)(l -x) + (m2 - (fc + p)2)x]"i+"2 

dxx^-^-^il- xY'^-^ d'^k 
[-fc2 -2xp-k- a;p2 + m2]"i+"2 ' 

dxx^-^-^il- x)"^-^ d'^k 



r(ni)r(n2) i7r''/2 J [m^ + x{l - x){-p'^) - A;'2]"i+"2 



we can integrate in k': 



r(ni +n2 " f) 



dx X 



no — l 



[l-x) 



rii — 1 



r(7ii)r(n2) 7o [m2 + a;(l -a;)(-p2)]"i+»2-d/2 
Now we shall use Mellin-Barnes representation 



1 



1 



(a + 5)" T{n) 27ri 



dz r(— z)r(?i + z) 



(5.16) 



Here the integration contour is chosen in such a way that all poles of r(- ■ ■ + z) (they are called 
left poles) are to the left of the contour, and all poles of r(- ■ ■ — z) (they are called right poles) 
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are to the right of it. It is easy to check (|5.16p : closing the contour to the right we get the 
expansion of the left-hand side in b/a; closing it to the left — the expansion in a/h. 
We continue our calculation: 



m 



d-2(ni+Ti2) I r + 



r(ni)r(n2) 27ri J_ 
m 



dzT{-z)T{ni+n2 + z) 



<i-2(«i+«2) 1 r+^oo r(-z)r(ni + z)r(n2 + z)r(?7,i + ?7,2 
' dz 



r(ni)r(n2) 2'Ki 



d I „ 



r(ni +n2 + 2z) 



-P 



This means that two massive lines can be replaced by one massless one (raised to the power 
—z) at the price of one extra integration in z: 



ni 




^d-2(«i+n2) 1 /■+'°°^^r(-z)r(7ii + z)r(?i2 + z)r(?ii + ?i2-rf/2 + z) 



r(ni)r(n2) 27ri 



r(ni +n2 + 2z) 



(5.17) 



This trick allows us to reduce this master integral to a single Mellin-Barncs integral. Using 
integration by parts, we can kill one of three lines in the left (integer) triangle, and calculate 
the integrand in F functions. This allows us to find several terms of its e expansion: 



r2(2£)r(3£- 1) 1 

4r(4e) 



-\-ioo 



dz 



r(l + z)r(l/2 + e + z)r(l + e + z)r{-2e - z)T{-e - z)T{^z) 



r(3/2 + £ + z)r(l-2e-z) 

TT^ 6C3 - StT 

9? ~ 



19 



' 11 ^4 10. , -2 
1 TT C-i H TT 

9e 270 3 ^ 9 



8, 82. 11 4 38, 65 2 



(5.18) 



This master integral has been evaluated in a closed form using Mellin-Barnes in a repre- 
sentation: 



r(l/2 - £)r(-£)r2(2e)r(l + e)T{3e - 1) 
' 4r(3/2 - e)r(4£) 



^( 2 +V^(l-£)-21og2 + 27B 



(5.19) 



This master integral can be written as a double Mellin-Barnes integral using (|5.17p . It 
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appears possible to calculate one integral: 
1 ^ V !■ 4^r(3/2 - e) 27ri 

+'°° ^_^ r(i + z)r (I - £ + z) r(g + z)r (^i + g - z) r (-§ + 2e-z) r(-z) 
-,oo ^ r(| + 2)r(e-z) 

r'^(l + ^^)^'^M-l + 2(41og2-7r-7)e + ---] . (5.20) 



5.2.4 Applications 

Feynman integrals considered here were used [20j for calculating the matching coefficients for 
the HQET heavy-quark field and the heavy-light quark current between the &-quark HQET with 
dynamic c-quark loops and without such loops (the later theory is the low-energy approximation 
for the former one at scales below TOc)- Another recent application — the effect of rric ^ on 
b ^ c plus lepton pair at three loops [22]. The method of regions was used; the purely soft 
region (loop momenta ^ m^) gives integrals of this type. Two extra terms of e expansion of the 
master integral of Sect. 15. 2?^ were required for this calculation which were not obtained in pp] . 
This was the initial motivation for |26j . 
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